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Abstract: Drought, which restricts the sustainable development of agriculture, ecological health, and 
social economy, is affected by a variety of factors. It is widely accepted that a single variable cannot fully 
reflect the characteristics of drought events. Studying precipitation, reference evapotranspiration (ET), 
and vegetation yield can derive information to help conserve water resources in grassland ecosystems in 
arid and semi-arid regions. In this study, the interactions of precipitation, ETo, and vegetation yield in 
Darhan Muminggan Joint Banner (DMJB), a desert steppe in Inner Mongolia Autonomous Region, China 
were explored using two-dimensional (2D) and three-dimensional (3D) joint distribution models. Three 
types of Copula functions were applied to quantitatively analyze the joint distribution probability of 
different combinations of precipitation, ET, and vegetation yield. For the precipitation—ETo dry—wet 
type, the 2D joint distribution probability with precipitationS245.69 mm/a or ET92959.20 mm/a in 
DMJB was approximately 0.60, while the joint distribution probability with precipitationS$245.69 mm/a 
and ETo=959.20 mm/a was approximately 0.20. Correspondingly, the joint return period that at least one 
of the two events (precipitation was dry or ET was wet) occurred was 2 a, and the co-occurrence return 
period that both events (precipitation was dry and ET» was wet) occurred was 5 a. Under this condition, 
the interval between dry and wet events would be short, the water supply and demand were unbalanced, 
and the water demand of vegetation would not be met. In addition, when precipitation remained stable 
and ET» increased, the 3D joint distribution probability that vegetation yield would decrease due to water 
shortage in the precipitation—ET» dry—wet years could reach up to 0.60—0.70. In future work, irrigation 
activities and water allocation criteria need to be implemented to increase vegetation yield and the safety 
of water resources in the desert steppe of Inner Mongolia. 
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1 Introduction 


Grassland is an important part of terrestrial ecosystems and plays an important role in regulating 
climate, conserving soil and water, improving soil quality, and maintaining biodiversity (Wang et 
al., 2022a). A desert steppe is the transition zone between grassland and desert, and it is the most 
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severe type of desertification among all grasslands. The desert steppe in Inner Mongolia 
Autonomous Region of China accounts for approximately 18.00% of the country's total grassland 
area (Li et al., 2013; Chu et al., 2014; Fang et al., 2016). The ecological environment of the desert 
steppe is fragile and extremely sensitive to precipitation changes (Zhan et al., 2007; Guo et al., 
2012; Song et al., 2017; Legesse et al., 2022). Strong evaporation and surface radiation aggravate 
the water sensitivity of the desert steppe vegetation (Bittelli et al., 2008; Jung et al., 2010; 
Armstrong et al., 2015). 

Precipitation and reference evapotranspiration (ETo) are two important, correlated random 
hydrological variables that characterize the relationship between regional water supply and plant 
demand (Li et al., 2017). At the same time, water is the main limiting factor affecting grassland 
ecosystems, which has been confirmed by controlled experiments and precipitation-gradient 
experiments (Bai et al., 2008; Yuan and Chen, 2009; Beier et al., 2012). In the desert steppe, 
monthly precipitation from May to August is decisive for annual precipitation. A study of the 
Inner Mongolia desert steppe reported that ETo was between 570.00 and 1674.00 mm/a for the 
period from 1961 to 2010, and it was highest in western Inner Mongolia (Wang et al., 2015). The 
aridity of the arid region is expected to increase in the 21‘ century (Andrés et al., 2017). Studying 
the response characteristics of vegetation yield (hereinafter abbreviated as yield) in the desert 
steppe to precipitation and ETo can not only reveal the quantitative risk of decreased yield 
because of water shortage under natural precipitation conditions, but also contribute to the 
understanding of ecological stability and yield security in the Inner Mongolia desert steppe. 

Many studies on arid and semi-arid ecosystems analyzed single precipitation event, or revealed 
the changes of a certain factor with respect to precipitation change (Yan et al., 2018; Yu et al., 
2019; Du et al., 2020). However, few researches have focused on the responses of yield to 
precipitation and ETo in the desert steppe. Moreover, it is also difficult to quantitatively analyze 
the multivariate distribution probability of these variables. Thus, to more accurately predict the 
influence of precipitation on yield, it is necessary to understand how yield responds to 
precipitation and ETo events. Copula functions are mathematical tools for solving multivariate 
probability problems and can reflect correlations between variables that are independent of the 
marginal distribution of random variables (Amirataee et al., 2020). Copula functions cannot limit 
the marginal distribution of variables. During the conversion process of the variables, the 
information of the original random variables will not change, especially when the variable 
dimension change is less than three. The multivariate probability model constructed from the 
Copula functions is simple and easy to apply (Xie et al., 2012). To date, Copula functions have 
been used extensively in hydrology-climate-related studies (e.g., Hao and Singh, 2015; Xu et al., 
2015; Zhang et al., 2015; Cai et al., 2021). Copula fitting was used to analyze the distribution, 
duration, and intensity of droughts (e.g., Amirataee et al., 2020). For example, Wu et al. (2015) 
used a Copula function to study the variations in hydrological drought frequency in southern 
China. Additionally, the drought risk can be qualitatively described and quantitatively evaluated 
through different combinations of precipitation and ETo, thereby improving irrigation planning 
and scientific deployment of water resources (e.g., Vo et al., 2020). 

To sum up, the three objectives of this study were as follows: (1) to obtain ETo using the 
Penman—Monteith formula, and choose and verify a suitable marginal distribution and a Copula 
function for precipitation, ETo, and yield of the desert steppe in Inner Mongolia; (2) to use a 
Copula function to quantitatively evaluate the two-dimensional (2D) and three-dimensional (3D) 
joint distribution probability of precipitation, ETo, and yield in the desert steppe; and (3) to 
explore the main factors that influenced the yield of the desert steppe, and provide measures to 
mitigate desertification and improve sustainable vegetation development in this region. 


2 Materials and methods 


2.1 Study area 
The study area is located in Darhan Muminggan Joint Banner (DMJB; 109°16'—-111°25’E, 
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41°20'—42°40'N), a typical desert steppe with an area of approximately 1.8110° km? in the 
central and western part of Inner Mongolia Autonomous Region, China. The area is a hilly desert 
steppe at the northern foot of Yinshan Mountain, and it is characterized by a temperate and 
semi-arid continental climate. The average annual precipitation is 253.45 mm, the annual average 
temperature is 4.12°C, and the average annual evaporation is 2480.57 mm (Song et al., 2017). The 
soil is dominated by brown soil and the vegetation is mainly perennial xerophytes. Species such 
as Artemisia frigida, Stipa krylovii, Stipa klemenzii, Stipa breviflora, and Cleistogenes songorica 
are the dominate plant species. The average vegetation coverage was 35.00% during 2002-2020. 
Data about the vegetation for this study were collected from the Yinshanbeilu Grassland 
Eco-hydrology National Observation and Research Station in DMJB (Fig. 1). 
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Fig. 1 Overview of land use/land cover in the study area (Darhan Muminggan Joint Banner; a), and photos 
showing Artemisia frigida (b) and Stipa krylovii (c) 


2.2 Data sources 


In this study, we collected daily meteorological data during 2002—2020 and annual data of yield in 
the main growing season (May—September) of 2002-2020. The daily precipitation, maximum 
temperature, minimum temperature, relative humidity, sunshine hours, solar radiation, water 
vapor pressure, and average wind speed data were provided by the China Meteorological 
Administration (http://cdc.cma.gov.cn). The data were collected and processed in accordance with 
the guide to the Global Observing System of the World Meteorological Organization and the 
technical specifications and standards for the weather observations of the China Meteorological 
Administration. The data of sunshine hours during July-December in 2020 were estimated from 
the average daily data of 2000-2019. The yield data of 2002-2020 were collected from an 
intensive pastoral region that covers an area of 1.33 km’. Artemisia frigida and Stipa krylovii were 
chosen as representative plant species to simulate yield in this study, as they dominated the 
composition of the plant community in the sampling plots and represented different functional 
groups (Zhao et al., 2006; Chen et al., 2013). In approximately mid-August every year (harvesting 
time), 5—10 plots (area of 1 mx1 m for each plot) with Artemisia frigida and Stipa krylovii as the 
dominant species were selected to survey, grid sample, measure, and evaluate for yield. One 
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sample was taken from one plot. From 2002 to 2020, a total of 167 grass samples were collected, 
and the dry weight of the aboveground part was measured. The data statistics and analysis were 
performed using R (Lucent Technologies, New Jersey, USA) and Excel software. The graphics 
were preprocessed by MATLAB R2018b software (MathWorks, Massachusetts, USA). The time 
ranges of the data for this study are shown in Table 1. 


Table 1 Description of the data used in this study 


Type Factor Unit Period 
Solar radiation kJ/(m?*-d) 2002-2020 
Minimum temperature °C 2002-2020 
Maximum temperature °C 2002-2020 
: Relative humidity % 2002-2020 
Meteorological data 
Water vapor pressure kPa 2002-2020 
Average wind speed m/s 2002-2020 
Precipitation mm/a 2002-2020 
Sunshine hours h/d 2000-2019 
Yield data Annual dry weight of yield per unit area g/(m°+a) 2002-2020 


Note: Yield, vegetation yield. 


2.3 Methods 


2.3.1 Calculation of ETo 


The Penman—Monteith formula was recommended by the United Nations Food and Agriculture 
Organization (FAO) in 1998 (Allen, 2000). It integrates data for elevation, latitude, wind speed, 
altitude, daily maximum temperature, daily minimum temperature, daily average temperature, 
daily average wind speed, daily average relative humidity, sunshine hours, and others. Then, ETo 
is calculated as follows (Zhang et al., 2012): 


EA E Oe 
T +273 i 


ay A +y(1+0.34u:) i 

where ETo is the reference evapotranspiration (mm/d); A is the slope of the saturated vapor 
pressure-temperature curve (kPa/°C); Rn is the net solar radiation of the plant canopy surface 
(MJ/(m*-d)); G is the soil heat flux (MJ/(m?-d)), and the daily G was calculated in this study 
(G=0); y is the hygrometer constant (kPa/°C); u2 is the wind speed at a height of 2 m (m/s); es and 
ea are the saturated water vapor pressure (kPa) and the actual water vapor pressure (kPa), 
respectively; and T is the average temperature (°C). 

2.3.2 Marginal distribution 


The Weibull (WEI), gamma (GAMMA), and Exponential (EXP) distributions were used to fit the 
marginal distributions of precipitation, ETo, and yield. The maximum likelihood method was 
applied to estimate the Copula function parameters, and then the Kolmogorov-Smirnov test (K-S 
test) was used to select the optimal marginal distribution (Song and Singh, 2010; Mohammadpour 
et al., 2012). 

The WEI distribution function is given by: 


Bey x? 
f(x) = (2) on( 2, (2) 


where f(x) is the marginal distribution function of x, where x is one of the variables (precipitation, 
ETo, or yield); a is the scale parameter; and b is the shape parameter. 
The GAMMA distribution function is given by: 


1 E 
f(x) re 'exo(-*), (3) 
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where T denotes the GAMMA function. 
The EXP distribution function is given by: 


f(x) = Aexp(-Ax), (4) 
where à is the scale parameter (A>0). 


2.3.3 Joint distribution model 

A Copula function does not limit the marginal distribution of variables. With the Copula model, 
several arbitrary marginal distributions can be connected to form a multivariate joint distribution 
probability model. Copula functions are based on the Sklar's theorem (Sklar, 1959). Random 
variables (Xi, X2,..., Xn (where n=3 in this study)) are continuous, and F1(X1), F2(X2), ..., Fn(Xn) 
(where F is the n-dimensional joint probability distribution function) are their marginal 
distribution functions. Then, there exists a Copula function (C) such that: 


F(X, Xz -o Xp) = CAQ, F(X), -> F, (X) (5) 
The K-S test statistic (D) that was used to select the optimal marginal distribution is defined as 
follows: 
k m—1 
D=maxs\e.- Je. - , (6) 
1<k<1 l l 


where cx is the Copula value of the joint observation sample (xz, yx) and mx is the number of joint 
observations that meet the conditions of x<x, and y<yx in the joint observation sample; k is the k™ 
number; and / is the total number (/=19 in this study). 

The Copula function is unique if the marginal distribution functions (F'1(x1), F2(x2), ..., FnQn)) 
are continuous. In this study, three commonly used functions, i.e., the Clayton Copula, 
Gumbel-Hougaard Copula, and Frank Copula functions (Lazoglou and Anagnostopoulou, 2019), 
were used to construct the 2D and 3D joint distributions of precipitation, ETo, and yield (Tables 2 
and 3). 


Table 2 Description of the two-dimensional (2D) Copula functions 


Copula function Function expression Condition 
1 
ae H(u,v)=(u?+yt-I) ¢ 20 
1 
Gumbel-Hougaard H(u,v) = exp He inu) + Cin vT ‘| 6>1 
1 -ou _4)( .- 4 
Frank H9)=—Le] io Me j 040 
et 


Note: The u and v are the marginal distribution functions of any two of precipitation, reference evapotranspiration (ETo), and yield; H(u, 
v) represents the joint distribution probability when both precipitation<uw (or ETo<w or yield<w) and precipitation<v (or ETo<v or 
yield<v) occur; and ĝis the Copula function parameter. 


Table 3 Description of the three-dimensional (3D) Copula functions 


Copula function Function expression Condition 
1 
Clayton H (yt ty) = (uy? +uy? +u;? -2) g 020 
2h 
Gumbel-Hougaard H(u,,Uy,U;) = exp{-[(-Inu,)? +(-Inu,)? + (-In us)’ | | o> 1 
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Note: The u1, u2, and u3 are the marginal distribution functions of precipitation, ETo, and yield, respectively; and H(u;, u2, u3) represents 
the joint distribution probability when precipitation<u), ETo<w2, and yield<w; all occur. 


Frank 
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The Copula function goodness-of-fit test is the premise for selecting the optimal Copula 
functions. Here, the fitting of the Copula functions was tested using the root mean square error 
(RMSE) and Akaike information criterion (AIC). The fitting test and the goodness-of-fit (Li et al., 
2013), with the minimum AIC as the best fit criterion, were calculated as follows: 


I g ; 
MSE = —— P.-P)’, 7 
vad! a =h) (7) 
RMSE =VMSE , (8) 
AIC = NIn(MSE)+2k , (9) 
Num(x, < x;, y; < y;)— 0.44 
P, (x; yi) = i aap g (10) 


where MSE is the mean square error; N is the total number of joint observations; Pe; is the 
empirical frequency; P; is the theoretical frequency; RMSE is the root mean square error; P.(x;, Vi) 
is the empirical frequency of the 2D joint distribution; Num(xj;<xi, yj<vi) indicates the number of 
joint observation value less than or equal to (x;, yi); and ź is the number of observation value. 

Num(x ; < Xi, Y; S Yi Z; S z;) -0.44 (11) 


t+0.12 
where P.(xi, yi, Zi) is the empirical frequency of the 3D joint distribution and Num(xj<xi, yi<yi, 
zj<zi) indicates the number of joint observation value less than or equal to (xi, vi, Zi). 
2.3.4 Joint distribution probability and return period 
Using the pf(wet) and pk(dry) (where pf and pk are the frequency values for dividing the annual 
precipitation or ETo into wet and dry, respectively) of 37.50% and 62.50%, respectively, as the 
frequency values for dividing the annual precipitation or ETo into wet, normal, and dry, we 
analyzed the interannual differences among them with the 2D Copula joint distribution model 
(Yan et al., 2007). We further determined the occurrences of wet, normal, and dry seasons, which 
can be divided into nine situations (Yan et al., 2007), as shown in Table 4. 


P, (Xis Yis Zi) = 


Table 4 Division of occurrences of wet, normal, and dry situations for precipitation-ETo 


Situation Frequency 
Precipitation—ET) wet—wet type P= PDX xy Y2V 5) 
Precipitation—ET,) wet—normal type Pr = P(X Z Xy Yy SY <p) 
Precipitation—ET» wet—dry type Pr = P(X 2X5, YSYp) 
Precipitation—ETy normal—wet type Pa = Pp <X < Xpy YZ Vy) 
Precipitation—ET, normal-normal type Ds = DX pp <X < Xpo Vor SY < py) 
Precipitation—ET) normal—dry type Do = P(X pp <X <X pps YS Vp) 
Precipitation—ET  dry—wet type Pp = P(X SX, YZY) 
Precipitation—ET» dry—normal type Ps = P(X SX p45 Y SY SY) 
Precipitation-ETo dry—dry type Py = P(X S Xp YS Yp) 


Note: ETo, reference evapotranspiration; pı—p9, the frequency of nine situations; X, the series of precipitation; Y, the series of (ET); pf and 
pk, the frequency values for dividing precipitation or ETo into wet and dry, respectively; x,¢and xpp, the values for dividing precipitation 
into wet and dry, respectively; y,r and y,,, the values for dividing ETo into wet and dry, respectively. 


Given that the main aim of this study was to evaluate the risk of decreased yield under the wet, 
normal, or dry condition of the Inner Mongolia desert steppe, we focused on the two joint 
distribution probabilities of the precipitation—ETo dry—wet type and their return periods: 

[o y(x, y)=P(X <xorY 2 y)=1- F; (y)+ F(x,y) 


i ; (12) 
Gy, yu y)= P(X < x,Y 2 y) = Fy (x)- Fy) 
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1 


Gy y(x, y) 
i ; (13) 


X,y 


xy 


Gy y(%y) 

where Gy, y(x, y) and P(X<x or Y>y) refer to the joint distribution probability when at least one of 
the two events that X does not exceed a certain value or Y exceeds a certain value occurs; G'y, y(x, 
y) and P(X<x, Y=y) refer to the joint distribution probability when two events that X does not 
exceed a certain value and Y exceeds a certain value occur; Fy(x) and Fy(y) are the marginal 
distributions of precipitation and ETo; F(x, y) is the joint probability distribution function of X<x 
and Y<y; Tx, is the joint return period; and 7‘, indicates the co-occurrence return period in 
which the two events occur simultaneously. 


3 Results 


3.1 Marginal distributions of precipitation, ETo, and yield, and division of precipitation 
and ET» 


The series of precipitation, ETo, and yield in DMJB from 2002 to 2020 are shown in Figure 2. 
ET» in DMJB was generally negatively correlated with precipitation and yield, while precipitation 
and yield were positively correlated with each other. The difference between precipitation and 
ETo fluctuated in the range of approximately 400.00—900.00 mm/a during the study period and 
was greatest in 2009, at 863.66 mm/a. 
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Fig. 2 Relationships between precipitation and reference evapotranspiration (ETo; a), precipitation and yield 
(vegetation yield; b), and ETo and yield (c) from 2002 to 2020 


With a significance level of a=0.05, precipitation, ETo, and yield in DMJB conformed to the 
GAMMA distribution (Table 5), which is a special form of Pearson type III distribution. When the 
significance level of the K-S test was 0.05 (n=19), the corresponding quantile was 0.30. The K-S 
test statistics that corresponded to precipitation, ETo, and yield for the GAMMA distribution were 
0.16, 0.11, and 0.17, respectively, which were all less than the corresponding quantile values, and 
the corresponding P values were all greater than 0.05. This indicates that the GAMMA 
distribution was appropriate and could provide values for dividing precipitation and ET» to wet, 
normal, and dry conditions in DMJB (Table 6). 


3.2 Joint distributions of precipitation, ET», and yield 


The maximum likelihood method was used to calculate the Copula functions parameters. The 
RMSE between the empirical frequency and theoretical frequency was calculated by the Copula 
functions, and the minimum AIC was used as the criterion for selecting the optimal fitting 
functions. The corresponding results are shown in Table 7. Because precipitation and yield were 
negatively correlated with ETo, only the Frank Copula function was suitable as their joint 
function. Precipitation and yield were positively correlated with each other. For the relationship 
between precipitation and yield, the Gumbel-Hougaard Copula function had the smallest AIC and 
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Table 5 Parameters of the univariate marginal distributions 


Characteristic Marginal distribution Parameter K-S test 
variable function Shape Scale Statistic P 
WEI 3.45 304.7600 0.20 0.38 
Precipitation GAMMA 13.18 0.0500 0.16 0.68 
EXP - 0.0036 0.45 0.04x107 
WEI 19.99 967.4000 0.13 0.85 
ETo GAMMA 288.14 0.3100 0.11 0.95 
EXP - 0.0011 0.59 8.59x107 
WEI 1.92 83.0500 0.19 0.45 
Yield GAMMA 4.11 0.0600 0.17 0.55 
EXP - 0.0100 0.36 0.01 


Note: WEI, Weibull; GAMMA, gamma; EXP, Exponential; K-S test, Kolmogorov-Smirnov test. K-S test statistic represents the 
maximum distance between the two distributions; the smaller the statistic, the smaller the gap between the two distributions and the 
more consistent the two distributions. The larger the P value, the more the null hypothesis cannot be rejected (P>0.05) and the more 
similarity between the two distributions. - indicates that there are no shape parameters. 


Table 6 Division values for dividing precipitation and ETo into wet and dry 


Wet (frequency of 37.50%) Dry (frequency of 62.50%) 
Precipitation (mm/a) ET  (mm/a) Precipitation (mm/a) ET» (mm/a) 
Division value 293.38 959.20 245.69 923.84 


Table 7 Parameter values and goodness-of-fit tests of the Copula functions 


Relationship of variables Copula function 0 AIC RMSE 
Precipitation-ETo Frank -2.71 —82.85 0.11 
Frank 4.13 62.95 0.18 
Precipitation-yield Clayton 1.03 -90.30 0.09 
Gumbel-Hougaard 1.84 -97.22 0.07 
ETo-yield Frank -1.11 —115.79 0.05 
Frank 0.09 -75.53 0.13 
Precipitation-ETo-yield Clayton 0.24 -72.85 0.14 
Gumbel-Hougaard 1.24 —69.86 0.15 


Note: 0, Copula function parameter; AIC, Akaike information criterion; RMSE, root mean square error. The Copula function with the 
smallest AIC and RMSE values was selected as the most appropriate probability distribution. 


RMSE values and the best goodness-of-fit. For ETo and yield, the Frank Copula function had the 
best goodness-of-fit, with AIC and RMSE values of -75.53 and 0.13, respectively. 

As shown in Figure 3, the empirical frequency and theoretical frequency of the joint 
distributions of precipitation, ETo, and yield in DMJB based on the Copula functions were 
uniformly distributed around the 45° diagonal line, indicating that the obtained Copula functions 
fitted well. The goodness-of-fit was best for precipitation—yield (Fig. 3c) and ETo-yield (Fig. 3d), 
with the R? values of 0.96 for both the empirical frequency and theoretical frequency, followed by 
precipitation—ET» (Fig. 3b), with the R? value of 0.88. However, the R? value for the relationship 
of precipitation—ETo—yield, at 0.76, was the smallest (Fig. 3a). These results show that the 
selected Copula functions were reasonable. 


3.3 Joint distribution probability of precipitation, ETo, and yield 


3.3.1 2D joint distribution probability 


The 2D joint distribution probability contours of precipitation—ETo, precipitation—yield, and 
ETo-yield, and the distributions of the actual values of precipitation, ETo, and yield during 
2002-2020 are plotted in Figure 4. The patterns in the trends of the 2D joint distribution 
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Fig. 3 Empirical frequency and theoretical frequency of joint distributions of precipitation—ETo-yield (a), 
precipitation—ETo (b), precipitation—yield (c), and ETo-yield (d) 
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Fig. 4 Two dimensional (2D) joint distribution contours of precipitation—ETo (a), precipitation—yield (b), and 
ETo-yield (c) in DMJB 


probability values for precipitation—ETo, precipitation—yield, and ETo-yield were the same. The 
actual years for precipitation—ETo joint distribution with precipitation<300.00 mm/a and 
ETo<1000.00 mm/a were 12 a in DMJB during 2002-2020, accounting for 63.16% of the total 
studied years, while the 2D joint distribution probability of precipitation—ETo with 
precipitation<300.00 mm/a and ET <1000.00 mm/a was 0.53 (Fig. 4a). The actual years for 
precipitation—yield joint distribution with precipitation<300.00 mm/a and yield<80.00 g/(m7-a) 
were 13 a, accounting for 68.42% of the total studied years, while the 2D joint distribution 
probability of precipitation—yield with precipitation<300.00 mm/a and yield<80.00 g/(m?-a) was 
approximately 0.53 (Fig. 4b). The actual years for precipitation—ETo joint distribution with 
ETo<1000.00 mm/a and yield<80.00 g/(m*-a) were 13 a in DMJB during 2002-2020, mainly 
concentrated in 2002-2016, and the 2D joint distribution probability of precipitation—ETo with 
ETo<1000.00 mm/a and yield<80.00 g/(m?-a) was approximately 0.53 (Fig. 4c). These results 
indicate that the ranges of the joint distributions of precipitation—ETo, precipitation—yield, and 
ETo-yield were relatively stable in DMJB during 2002-2020. The probability that yield would be 
<80.00 g/(m?-a) was approximately 0.50. 
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3.3.2 Frequency analysis of the 2D joint distributions of the wet, normal, and dry conditions for 
precipitation—ETo 

For the nine combinations of wet, normal, and dry conditions, the combination of the dry-wet 
type for precipitation—ETo was the most unfavorable for irrigation scheduling. Figure 5 shows the 
2D joint distribution contours of precipitation—ETp and the distributions of the actual precipitation 
and ETo values in DMJB during 2002-2020. Specifically, Figure 5a indicates the joint distribution 
probability that at least one of the two events occurred, i.e., precipitation was dry or ETo was wet, 
and Figure 5b shows the joint distribution probability that both events occurred, i.e., precipitation 
was dry and ETo was wet. Overall, the patterns of changes in the joint probability values were the 
same, and they all increased with the increase of precipitation and the decrease of ETo. The joint 
distribution probability for the precipitation-ETo dry—wet type with the actual value of 
precipitation<245.69 mm/a or ETo=959.20 mm/a was approximately 0.60 during 2002-2020, and 
there were 12 a in this range, accounting for 63.16% of the total studied years. However, the joint 
distribution probability for the precipitation—ETo dry—wet type with precipitation<245.69 mm/a 
and ETo2959.20 mm/a was approximately 0.20, and there were 3 a in this range, accounting for 
15.79% of the total studied years. 

The joint return period that at least one of the two events (precipitation was dry or ETo was 
wet) occurred was 2 a, and the co-occurrence return period that both events (precipitation was dry 
and ETo was wet) occurred was 5 a. These values were calculated from Equation 11 and were 
based on the joint distribution probability of the precipitation—ETo dry—wet type in DMJB. The 
water supply and demand imbalance of DMJB is more likely when there is natural precipitation. 
Under this condition, the interval between dry and wet events would be short, and the risk of 
vegetation water demand not being met is greater. 


(a) Precipitation was dry or ET, was wet (b) Precipitation was dry and ET, was wet 
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Fig. 5 2D joint distribution probability contours of precipitation—ETo with at least one of the two events 
(precipitation was dry or ETo was wet) occurred (a) and with both events (precipitation was dry and ETo was wet) 
occurred (b) 


3.3.3 3D joint distribution probability 

The analysis in Table 7 indicated that the Frank Copula function could be applied to the 3D joint 
distribution of precipitation, ETo, and yield when precipitation<u; (where wu; is the marginal 
distribution function of precipitation), ETo<w2 (where u2 is the marginal distribution function of 
ET»), and yield<u3 (where u3 is the marginal distribution function of yield) all occurred in DMJB. 
The 3D joint distribution probability of precipitation, ETo, and yield is shown in Figure 6. The 
average yield per unit area was 73.05 g/(m?.a) in DMJB during 2002-2020. When precipitation 
was dry and ET» was wet, as in 2004, 2005, and 2009, the average yield per unit area was at a low 
level (48.74 g/(m?.a)), and the 3D joint distribution probability of precipitation-ETo—yield with 
precipitation<245.69 mm/a, ET <959.20 mm/a, and yield<48.74g/(m?.a) was 0.07. When 
precipitation was dry or ETo was wet, the average yield per unit area was 68.24 g/(m?.a), and the 


naXiv:202304.00936v1 


i 


ch 


ChinX ivê /ERAF! 


LI Hongfang et al.: Responses of vegetation yield to precipitation and reference... 487 


3D joint distribution probability of precipitation-ETo-yield with precipitation<245.69 mm/a, 
ETo<959.20 mm/a, and yield<68.24 g/(m°-a) was 0.12. Of these, when precipitation was stable 
and ETo increased, the 3D joint distribution probability that yield would decrease due to water 
shortage in a year with the precipitation—ET» dry—wet type could reach 0.60—0.70. In addition, as 
shown in Figure 6, when ETo was lower than 900.00 mm/a, the 3D joint distribution probability 
that yield would decrease was less than 0.30, indicating that the yield of the desert steppe was 
sensitive to precipitation and other factors such as temperature, wind speed, and sunshine hours. 


3D joint distribution probability 
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Fig. 6 Three dimensional (3D) joint distribution probability of precipitation, ETo, and yield in DMJB 


4 Discussion 


During the crop growth period, precipitation and ET» are the main factors that cause agricultural 
drought. Batsukh et al. (2021) used ET» to predict the evapotranspiration of a specific biological 
community in the desert steppe of Inner Mongolia in the absence of meteorological data. Han et 
al. (2015) studied the temporal and spatial dynamics of climate and vegetation phenology in the 
desert steppe of Inner Mongolia and found that the primary productivity decreased from southeast 
to northwest along the precipitation gradient. Precipitation can affect vegetation growth when it is 
transformed into soil water through interception by the vegetation canopy (Sugita et al., 2015). 
However, the previous studies are relatively one dimensional and lack of multivariate research on 
the desert steppe (Zhu et al., 2022). The joint distribution probability for the precipitation—ETo 
dry—wet type with the actual value of precipitation<245.69 mm/a or ETo=959.20 mm/a was 
approximately 0.60 during 2002—2020 in our study. 

On the basis of previous research, we quantitatively analyzed and predicted the 2D and 3D 
joint risk probabilities of precipitation, ETo, and yield of the desert steppe in Inner Mongolia. The 
results show that for the precipitation—ETo dry—wet type, the 2D joint distribution probability with 
the actual precipitation<245.69 mm/a or ETo=959.20 mm/a in DMJB was approximately 0.60, 
while the 2D joint distribution probability with precipitation<245.69 mm/a and ETo=959.20 mm/a 
was approximately 0.20. The joint return period and the co-occurrence return period were 
calculated as 2 and 5 a, respectively (Fig. 4), indicating that drought risk is higher in the 
precipitation—ETp dry—wet years, which could increase the degradation of grasslands in Inner 
Mongolia (Hu et al., 2019). This result is consistent with the finding that evapotranspiration 
would further increase under the background of global warming (Li et al., 2016). The desert 
steppe was also affected by changes in ETo (Figs. 3c and 6), which is consistent with the fact that 
ETo in the temperate desert steppe of Inner Mongolia can be well estimated using only 
meteorological variables and the Penman—Monteith equation (Yang and Zhou, 2011). Previous 
studies also demonstrated that the effect of small precipitation events (0.00-5.00 mm) on plant 
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assimilation was limited (e.g., Song et al., 2017). Our results showed that when ETo was lower 
than 900.00 mm/a, the 3D joint distribution probability of decreased yield would be less than 0.30 
(Fig. 6), suggesting that the yield of the desert steppe was sensitive to precipitation and other 
factors such as temperature, wind speed, and sunshine hours (Wang et al., 2022b). Therefore, 
artificial precipitation, reduction of grazing, and other measures are needed in years with the 
precipitation—ETo dry—wet type to ensure the sustainable development of the desert steppe, and 
adaptability to drought is also essential for the growth of grassland species. 

This study plays an important role in analyzing the probability of dry and wet events during the 
crop growth periods in the desert steppe of Inner Mongolia and guiding agricultural drought 
control. However, there were some limitations in data processing in our work. The vegetation data 
were limited. As Artemisia frigida and Stipa krylovii are the dominant plant species in DMJB, rare 
varieties were ignored in vegetation collection. In addition, although the suitable Copula functions 
selected by the marginal distributions of precipitation, ETo, and yield had a good fit (Fig. 3; Table 
7), there were certain deviations in the software precision. In the future, related studies on the 
desert steppe should consider other factors, such as biodiversity. It is also worth exploring the 
feedback mechanisms between soil organisms and hydrological cycle in arid and semi-arid 
regions. 


5 Conclusions 


This study focused on the response of yield to precipitation and ETo in wet and dry years in a 
typical desert steppe (DMJB) in Inner Mongolia. We also examined how precipitation, ETo, and 
yield interacted in DMJB during 2002-2020. The three variables all obeyed the GAMMA 
distribution, and the joint distribution models of precipitation—ETo—yield, precipitation—ETo, and 
ETo-yield established by the Frank Copula function fitted well. The Gumbel-Hougaard Copula 
function exhibited the best fit for the joint distribution model of precipitation—yield. The R? values 
of the empirical frequency and theoretical frequency were all higher than 0.70, and the 3D joint 
distribution probability of precipitation—ETp—yield with precipitation<245.69 mm/a, ETo<959.20 
mm/a, and yield<48.74g/(m7-a) was 0.07. The joint distribution probability values of precipitation 
and ET» all increased with the increase of precipitation and the decrease of ETo in DMJB. Under 
natural precipitation, there was a relatively high risk in DMJB that the water supply would not 
meet the water demand of the vegetation, and the interval time between dry and wet events was 
short. For the precipitation—ETo dry—wet type, the 2D joint distribution probability with 
precipitation<245.69 mm/a or ETo=959.20 mm/a in DMJB was approximately 0.60, while the 
joint distribution probability with precipitation<245.69 mm/a and ET 9=959.20 mm/a was 
approximately 0.20. Further, the joint return period and the co-occurrence return period were 2 
and 5 a, respectively. Irrigation activities and water allocation criteria need to be implemented to 
ensure the success of agriculture and animal husbandry in the precipitation—ETo dry—wet years. 

In the future, the research scope can be expanded to the whole Inner Mongolia grassland, then 
the location-varied patterns can be deduced through the spatial variations, and the joint 
distributions of precipitation, ETo, and yield can be analyzed at different locations. 
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